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We consider advection of a passive scalar 9(t, r) by an incompressible large- 
scale turbulent flow. In the framework of the Kraichnan model the whole PDF's 
(probability distribution functions) for the single-point statistics of 9 and for 
the passive scalar difference 9{r{) — 9(r?) (for separations t*i — r2 lying in the 
convective interval) are found. 



INTRODUCTION 

> 

Q\ \ We treat advection of a passive scalar field 9(t, r) by an incompressible turbulent flow, the role of 
the scalar can be played by temperature or by pollutants density. The velocity field is assumed to 
contain motions from some interval of scales restricted from below by L v . A steady situation with a 
permanent random supply of the passive scalar is considered. We wish to establish statistics of the 
passive scalar 9 for scales that are less than both the scale L v and the pumping scale L, and larger 
than the diffusion scale fdif (for definiteness we assume that L < L v ). Such convective interval of 
£>y scales exists if the Peclet number Pe = L/r&t is large enough, we will assume the condition. Since 
all scales from the convective interval are assumed to be smaller than L v we will say about advection 



by a large-scale turbulent flow. The problem is of physical interest for the dimensionalities d = 2,3, 



but formally it can be treated for an arbitrary dimensionality d of space. Below we will treat d as a 
parameter. Particularly, all the expressions will be true for a space of large dimensionality d. 
! Description of a small-scale statistics of a passive scalar advected by a large-scale solenoidal velocity 
field is a special problem in turbulence theory. This problem was treated consistently from the very 
beginning and some rigorous results have been obtained which is quite unusual for a turbulence 
problem. Batchelor |l] examined the case of external velocity field being so slow that it does not 
change during the time of the spectral transfer of the scalar from the external scale to the diffusion 
scale. Then Kraichnan || obtained plenty of results in the opposite limit of a velocity field delta- 
correlated in time. The pair correlation function of the passive scalar (9(r)9(0)) was found to 
be proportional to the logarithm ln(L/r) and the pair correlation function of the passive scalar 
difference ([9(r) — 9(0)] 2 ) was found to be proportional to ln(r/r^) in both cases. The assertions 
are really correct for any temporal statistics of the velocity field Thus we are dealing with 

the logarithmic case which is substantially simpler than cases with power-like correlation functions 
usually encountered in turbulence problems 0-0]. 

Now about high-order correlation functions of the passive scalar. As long as all distances between 
the points are much less than L, the 2n-point correlation functions of 9 are given by their reducible 
parts (that is are expressed via products of the pair correlation function) until n ~ ln(L/r) where r 
is either the smallest distance between the points or r^f depending on what is larger M . The reason 
for such Wick decoupling is simply the fact that reducible parts contain more logarithmic factors 
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(which are considered as large ones) than non-reducible parts do. Consistent calculations of the 
fourth-order correlation function of the passive scalar at d — 2 || confirm the assertion. Therefore 
e.g. the single-point PDF of 9 has a Gaussian core (that describes the first moments with n < In Pe) 
and a non-Gaussian tail (that describes moments with n > In Pe). The tail appears to be exponential 
(see [0,|J). The same is true for the passive scalar difference A9 = 9(r) — 9(0) where instead of 
InPe we should take ln(r/rdif). The tails do not depend on InPe or on ln(r/rdif) and contain only 
coefficients depending on the statistics of the advecting velocity. 

Correlation functions of the passive scalar can be written as averages of integrals of the pumping 
along Lagrangian trajectories (see e.g. [[|). For example, the pair correlation function (9(r)9(0)) 
is proportional to an average time needed for two points moving along Lagrangian trajectories to 
run from the distance r to the distance L. Generally, correlation functions of the passive scalar are 
determined by the spectral transfer that is by an evolution of Lagrangian separations up to the scale 
L. For the large-scale velocity field Lagrangian dynamics is determined by the stretching matrix 
c"a/3 = V/jfc and, consequently, the statistics of the matrix determines correlation functions of the 
passive scalar. For example, the coefficient at the logarithm in the pair correlation function of the 
passive scalar is -P2/A EH1 where P 2 is the pumping rate of 9 2 and A is Lyapunov exponent that is 
the average of the largest eigen value of the matrix a. The coefficients in the exponential tails are 
more sensitive to the statistics of a, particurlaly they depend on the dimensionless parameter At Q 
where r is the correlation time of a. The motion of the fluid particles in the random velocity field 
resembles in some respects random walks, but one should remember that correlation lengths of both 
the advecting velocity and of the pumping are much larger than scales from the convective interval 
we are interested in. Thus the situation is opposite to one usually occurring in solid state physics, 
where e.g. random potential is shortly correlated in space. 

Since really ln(L/r) is not very large it is of interest to find the whole PDF's for the single-point 
statistics of 9 and for the passive scalar difference A9. It is possible to do for the Kraichnan short- 
correlated case Ar C 1 when the statistics of a can be regarded to be Gaussian. An attempt 



to do this was made in |T0| , |TT| in terms of the statistics of the main eigen value of the matrix a. 
Unfortunately the scheme works only for the dimensionality d = 2 where the matrix a has a single 
eigen value. This fact was noted in the work [T2j where also the correct coefficient in the exponential 



tails for an arbitrary dimensionality of space d was found. Here, we develop a scheme enabling to 
obtain the whole PDF's for arbitrary d. The scheme is also interesting from the methodical point of 
view. For example, its modification enables one to calculate the statistics of local dissipation ||13|| . 

The paper is organized as follows. In Section | we find a path integral representation for the 
simultaneous statistics of the passive scalar. In Section |T| we analyze the generating functional 
for correlation functions of the passive scalar in the convective interval of scales. Using different 
approaches we obtain the functional and establish the applicability conditions of our consideration. 
In Section |H| we find explicit expressions for the single-point PDF and for the PDF of the passive 
scalar difference. In Conclusion we shortly discuss the obtained results. 

I. GENERAL RELATIONS 

The dynamics of the passive scalar 9 advected by the velocity field v is described by the equation 

d t 9 + vW9 - kV 2 9 = <f). (1.1) 

Here, the term with the velocity v describes the advection of the passive scalar, the next term is 
diffusive (k, is the diffusion coefficient) and <fi describes a pumping source of the passive scalar. Both 
v(t, r) and 4>(t, r) are supposed to be random functions of t and r. We regard the statistics of the 
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velocity and of the source to be independent. Therefore, all correlation functions of 9 are to be 
treated as averages over both statistics. 



A. Simultaneous Statistics 



The source is believed to possess a Gaussian statistics and to be 5-correlated in time. The 
statistics is entirely characterized by the pair correlation function 

(<f>{h, n)<f>{t 2 , r 2 )> = 8(h - t 2 ) x (In - n|) , (1.2) 

where we assumed that the pumping is isotropic. The function x( r ) * s accepted to have a charac- 
teristic scale L which is the pumping length. We will be interested in the statistics of the passive 
scalar on scales much smaller than L. 

Simultaneous correlation functions of the passive scalar 9 can be represented as coefficients of the 
expansion over y of the generating functional 

Z(y) = (expjky J dr 0(r)8(O, r)}^ , (1.3) 

where (3 is a function of coordinates and angular brackets designate averaging over both the statistics 
of the pumping and over the statistics of the velocity v. The generating functional Z(y) contains 
the entire information about the simultaneous statistics of the passive scalar Particularly, knowing 
Z(y) one can reconstruct the simultaneous PDF of the passive scalar, the problem is discussed in 
Section |TJ. 

If characteristic scales of (3 in (|1.7|) are much larger than the diffusion scale r^f then it is possible to 



neglect diffusivity at treating the generating functional (|1.3|) . Then the left-hand side of the equation 
( |1.1|) describes the simple advection and it is reasonable to consider a solution of the equation in 
terms of Lagrangian trajectories g(t) introduced by the equation 

d t Q = v(t, g) . (1.4) 

We will label the trajectories by r, those are positions of Lagrange particles at t — 0: g(0, r) = r. 
Next, introducing 9(t, r) = 9(t, g) we rewrite the equation (|1.1|) as d t 9 = 0, that leads us to 

0(0, r)= f° dt<f>(t, ). (1.5) 

J — oo 

Here we have taken into account that at t = the functions 9 and 9 coincide. Starting from ( |1 . 5[ ) 
and using Gaussianity of the pumping statistics we can average explicitly the generating functional 
( |Op over the statistics. The result is 



Z(y) = ( exp 



2 



dtU 



;i.6) 



U = J dr 1 dr 2 /5(r 1 )/3(r 2 )x(|^ 1 - Q 2 \) , (1.7) 

where angular brackets mean averaging over the statistics of the velocity field only. 

Being interested in the single-point statistics of 9 we should take (3(r) = 5(r). But it is impossible 
since we have neglected diffusivity. We will take (3{r) = 5\(r) instead where the function 5a(v) 
tends to zero at Ar > 1 fast enough and is normalized by the condition 
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dr 5 A (r) = 1 . 

Then the generating functional (1.6) will describe the statistics of an object 

6 A = f dr5 A (r)6(r), (U 



smeared over a spot of the size A -1 . If rdifA <C 1 then the statistics of the object is not sensitive 
to diffusivity. From the other hand if AL ^> 1, then knowing correlation functions of 9 A , we can 
reconstruct single-point statistics due to the logarithmic character of the correlation functions. To 
obtain single-point correlation functions one should substitute simply A — > r~j^ into the correlation 
functions of 6 A . The above inequalities Ar^f -C 1 and AL ^> 1 are compatible because of Pe ^> 1. 
If we are interested in the statistics of the passive scalar differences in points with a separation r 
(where r 3> raif) then instead of 5 A (r) we should take 

P(r) = 5 A (r - t-q/2) - 5 A (r + r /2) . (1.9) 

Then the generating functional ( |1.6| ) will describe the statistics of an object 

A6 A = 9 A (r /2)-e A (-r /2). (1.10) 

Again, correlation functions of the passive scalar differences can be found from correlation functions 
of A9 A after the substitution A — > r^. 



B. Path Integral 

Below, we treat advection of the passive scalar by a large-scale velocity field, that is we assume 
that the velocity correlation length L v is larger than the scales from the convective interval. Then 
for the scales one can expand the difference 

v a (ri) - v a (r 2 ) = a a/3 (t)(r ll3 - r 2 p) , a af3 = V ' pv a . (1.11) 

Here a a p{t) can be treated as an r-independent matrix field. Then the equation (|L4|) leads to 

9 t {gi,a ~ Q2,a) = CTap{t){Ql,p ~ Q2,p) ■ (1-12) 

A formal solution of the equation ( |1.12j ) is 

Ql,a - Q2,a = W a /3(ri,p ~ r 2 fi) , 

d t W = aW, W = Texp (- [ dta] , (1.13) 



where T means an antichronological ordering. Note that Det W — 1, those property is a consequence 
of Tro" = and the initial condition W = 1 at t = 0. The Lagrangian difference in ( |1.7| ) is now 
rewritten as 



If?x - 2 | = V( r i« - r 2a )B aP {r ip - r w ) , B = W T W , (1.14) 

where the subscript T designates a transposed matrix. Note that Det-B = 1 since Det W = 1. 

The generating functional Z(y) (|1.6| ) can be explicitly calculated in the Kraichnan case when 
the statistics of the velocity is 5-correlated in time. Then the velocity statistics is Gaussian and is 
entirely determined by the pair correlation function which in the convective interval is written as 
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(v a (ti, ri)vp(t2, r 2 )) = 5(ti - t 2 ) [V 5 a/ 3 - K a p{r\ - r 2 )] , 
K, a p{r) = D(r 2 5 al3 - r a rp) + — ^ D 8 af3 r 2 . 



;i.i5) 

1.16) 



Here Vo is a huge r-independent constant and D is a parameter characterizing the amplitude of the 
strain fluctuations. The structure of (|1.16|) is determined by the isotropy and space homogeneouty 
accepted, and by the incompressibility condition div v = 0. Then the statistics of a is Gaussian and 
is determined by the pair correlation function which can be found from ( p.. 15 , 1.16 ): 



{cr af 3{ti)a^ v {t 2 )) = D[(d+ l)5 atl 5p v - 5 av 8p^ - 8 a p8^ v ] 6(ti - t 2 ) 



;i.i7) 



Note that the correlation function ( |1 . 1 7| ) is r-independent, as it should be. We see from ( |1.17| ) that 
the parameter D characterizes amplitudes of a fluctuations. 

Averaging over the statistics of a can be substituted by the path integral over unimodular matrices 
W(t) with a weight exp (il). The effective action I = J dt Co is determined by ( |1 . 1 7|) : 



i£ 



2d(d + 2)D 



(d + l)Tr(a T a) + Tra 5 



:i.i8) 



Then the generating functional (|1.7|) can be rewritten as the following functional integral over uni- 
modular matrices 



Z(y) 



VW exp 



dt 



u 



U = J &ridr 2 f3(ri)f3(r 2 )x \j (r la - r 2a )B af) {r ip - r 2/3 ) 



;i.i9) 

1.20) 



Here, we should substitute a = d t W(W) 1 and remember about the boundary condition W = 1 at 
t = 0. 

Some words about the 'potential' U (|1.7|) figuring in ( |1.20|) . The characteristic value of t*i — r 2 in 
the integral ( |1.7| ) is of order A -1 for (3(r) = 5\(r). Since we assume AL ^> 1 then for a single-point 
statistics U ~ P 2 , where P 2 = x(0), if B is not very large. Particularly it is correct at moderate 
times \t\ since B = 1 at t = 0. With increasing \t\ the argument of x i n (|1-20|) grows and U 
tends to zero when the argument of x becomes greater than L. For the passive scalar difference 
when (3 is determined by ( |1.9|) the situation is a bit more complicated. Then U is a difference 
of two contributions. The first contribution behaves like for a single-point statistics. The second 
contribution contains x with the argument determined by r% — r 2 w ±r . Then at t = the meaning 
of the second contribution is determined again by P 2 but it turns to zero at increasing \t\ earlier 
than the first contribution. 

Thus we reduced our problem to the path integral ( |1.19| ) that is to the quantum mechanics with 
d 2 — 1 degrees of freedom. Nevertheless to solve the problem we should perform an additional 
reduction of the degrees of freedom. The conventional way to do this is passing to eigen values, say, 
of the matrix B figuring in (|1.20| ) (see e.g. [[14]]) and excluding angular degrees of freedom. Just 
this way was used by Bernard, Gawedzki and Kupiainen in [0. Then the authors using known 

E3 



have found the 



facts about the quantum mechanics associated with the eigenvalues (see e.g. 
coefficient in the exponential tail in the single-point PDF of 9. Unfortunately this way is not very 
convenient to find the whole PDF. To do this we will use a special representation of the matrix W in 
the spirit of the nonlinear substitution introduced by Kolokolov |Tj|. That is a subject of the next 
subsection. 
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C. Choice of Parametrization 



To examine the generating functional Z{y) we will use a mixed rotational-triangle parametrization 

W = RT, B = f T f 1 (1.21) 



where R is an orthogonal matrix and T is a triangular matrix: Ty = at % > j. The parametrization 
( |1.21| ) is the direct generalization of the 2d substitution suggested in [17]. Note that DetT = 1 since 
Det W = 1. Note also that the matrix B introduced by (|1.14| ) does not depend on R as is seen 
from ( p.. 21 ). That is a motivation to exclude the matrix R from the consideration performing the 
integration over the corresponding degrees of freedom in the path integral ( |1.19|) . A jacobian appears 
at the integration. To avoid an explicit calculation of the jacobian (which needs a discretization over 
time and then an analysis of an infinite matrix [|l(J) we will use an alternative procedure described 
below. 

Let us examine the dynamics of the matrix T. It is determined by the equation 



i<k<j 

following from ( |1.13| , |L~2~T|) . Here we used the designation 

E = R T aR . 

Next introducing the quantities 

Tu = exp(pi) . 
we rewrite the equation ( |1.22| ) as 

dtPi = , 

d t rj i:j = (Ey + Eji) exp(Pj 



Tij = exp(pi)^ , if i < j 



Pi 



+ eX P(Pfc - Pi)Vkj 

i<k<j 



Comparing (|1.13|) with ( |1.21| ) one can find the following expression for A = R T dtR 



A 



ij 



E {j if i > j 



A; 



-Eji if i < j 



(1.22) 

(1.23) 

(1.24) 

(1.25) 
(1.26) 

(1.27) 



One can easily check that the irreducible pair correlation function of has the same form as for 
a"jj that is determined by ( |1.17|) : 



(Ey (ti)£ mn (t 2 )} = D[(d + l)5 im 5 jn - 5 in 5 jm - <%<5 mn ]#(*i - h) ■ 
Besides, the average value of Ey is nonzero ||10|| : 



i 3 ■ 



1.28) 



1.29) 



Nonzero averages of Ey are related to Lyapunov exponents (not only the first one) [|ll| (for our 
model see also [0). To obtain ( |1 . 29| ) one should first solve the equation A = R T d t R for R on a 
small interval r: 



R(t + t) w R(t) 



1 + 



t+T 



dt'A{t') 



G 



Then with the same accuracy we get from the expression ( |1.23j ) 



£(t + r) ss R T (t)a(t + T)R(t) 



t(t + r) 



t+T 



dt 1 A{t') 



1.30) 



The average value of S arises from the second term in the right-hand side of ( |1.30|) . The explicit 
form of the average can be found using 

/ f t+T i \ D 

(Y,ij(t + r) dtE mn (t)) = —[(d + l)S im S jn - 5 in 5 jm - 5ij5 mn ] . (1.31) 
Here we utilized Eq. ( |1.28|) and substituted the integral 

rt+r 

J dt' 8{t + T- t') 

by 1/2. The reason is that really the correlation function of a has a finite correlation time and 
therefore S(t) (representing this correlation function) should be substituted by a narrow function 
symmetric under t — > —t. Then we will get 1/2. Expressing in ( |1.30| ) A via £ from ( |1.27| ) and 
calculating its average using ( |1.31| ) we get the answer ( |1.29| ). 

The expressions ( |1 . 25| . p~726| . |l . 28|JIT29| ) entirely determine the stochastic dynamics of and rjij. 
Using the conventional approach |2T)|-pM| correlation functions of these degrees of freedom can be 
described in terms of a path integral over p i; rjij and over auxiliary fields which we will designate as 
rrii and (i < n). This integral should be taken with the weight exp(i / dt C) where the Lagrangian 
is 

/ \ 21 



a=l 



dtPa + D 



d(d-2a+ 1) 



+ 



LD 



+iDdJ2 ex P( 2 Pj ~ 2 Pi)l j2 ij + 2iDd PijPikexp(2p k ~2pi)r] kj 

i<j i<k<j 

+ PijdtVij + lD d PimPinVkmVkn exp(2p fc - 2pi) . 

i<j i<k<m,n 



1.32) 



Since the matrix B in accordance with Q1.21 ) does not depend on R it is enough to know the statistics 
of p a and r)ij to determine the average (|1.6|). Therefore, instead of (|1.19|) we get 

,2 



Z(y) = J V pVr/VmV p, exp 



dt 



(1.33) 



Here U is determined by ( |1.20| ) where (|1.21 , L2i ) should be substituted. 

Thus we obtained the expression for the generating functional ( |1.3|) in terms of the functional 
(path) integral which is convenient for an analysis which is presented in the subsequent section. 



II. GENERATING FUNCTIONAL 

Here, we are going to calculate the generating functional (|1.3| ) for a single-point statistics of 9 that 
is of the object (|L8|) corresponding to (3{r) = 8\(r) and also the statistics of the difference that 
is of the object (|1.10|) corresponding to (jl.9Q . The starting point for the subsequent consideration 
is the expression (|1.33j ). There are different ways to examine Z(y). We will describe two schemes 
leading to the same answer but carrying in some sense a complimentary information. We believe 
also that the consideration of the different schemes is useful from the methodical point of view. A 
modification of the second scheme is presented in Appendix. 
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A. Saddle-Point Approach 



The first way to obtain the answer for the generating functional ( |1.3| ) is in using the saddle-point 
approximation for the path integral (|1.33| ). The inequalities justifying the approximation are AL ^> 1 
for the object ( |1.8| ) and Ar ^> 1 for the object (|1.10| ). 

As we will see large values of the differences pi — pk (i < k) will be relevant for us. At the condition 
fluctuations of r\ and p are suppressed and it is possible to neglect the fluctuations. Therefore we can 
omit the integration over r\ and p in ( |1.33|) substituting r\ = p = into (|1.32|) . After that we obtain 
a reduced Lagrangian. Introducing an auxiliary field s we can rewrite the reduced Lagrangian as 



i£ r 



d 

i V] m 

a=l 



OtPa + D 



d(d+l- 2a) 



+ s 



Dd x , 2 

> m n 

2 V 



s 

2D 



(2.1) 



Now, to obtain Z(y) one should integrate the exponent in ( 1.33| ) (with C r 
examine Q2.ip it is convenient to pass to new variables 4> a = O a bpb and fh a 
orthogonal matrix. We make the following transformation 



d(d 2 - 1) 



[{d - l) Pl + {d- 3)p 2 + . . . + (1 - d) Pa 



over p a , m a and s. To 
= O a bmb where O is an 



(2.2) 



Then the expression (2.1) will be rewritten as 



LC r = i m a dt4> a + iVa: 



m d s 



s 

2D 



Dd 




mi . 



(2.3) 



The Lagrangian ( |2.3|) is a sum over different degrees of freedom. The dynamics of 0i is ballistic 
whereas the dynamics of <p a for a, > 1 is purely diffusive. We will see that times determining the 
main contribution to the generating functional are large enough so that for the relevant region 
0i 3> 4>a- Therefore the potential U ( |1.20| ) depends practically only on <pi and it is possible to 
integrate explicitly over s, <p a and rh a for a > 1. After that we stay only with one degree of freedom 
which is described by the Lagrangian 



LCi = irfii 




Dd 



T-,2 



-m 



(2.4) 



Neglecting all (fi a for a > 1 and performing the transformation inverse to (|2.2j) we obtain 



Pi 



3(d 



\ d(d+l) 



Pa 



d-2a + l 
d- 



-Pi 



(2.5) 



We will soon see that the characteristic value <pi ^> 1. Therefore the characteristic value of e pl is 
much larger than other e pa and we conclude that the potential U depends really only on p\. For 
the case of the single-point statistics the characteristic value of the difference r\ — r 2 in (|1.20|) is 
A -1 . Then it follows from ( |1.21| , |L2"4| ) that the potential U falls down from P 2 to zero near the point 
pi = ln(LA) that is near the point <fti = <p\ where 



d(d+l) 
\ 3(d- 1) 



In (LA) 



(2.6) 
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For the difference the potential increases from zero to 2P 2 at (j)\ = 4>b, where 



d(d + 1) 
\ 3(d- 1) 



ln(L/r ) 



(2.7) 



and then falls down from P 2 to zero near the point (pi = <p\. The expressions (|2.6|J2.7| ) determine 
the characteristic values of <pi which are really large since LA > 1 or L/r > 1 what justifies our 
conclusions. 

Now we can employ the saddle-point approximation: 



In Z{y) 



dt 



(2.8) 



inst 



where we should substitute solutions of the extrema conditions which we will call instantonic equa- 
tions. The instantonic equations which can be found as extrema conditions for LCi — y 2 U/2 are 



Dd did 2 



1} 



-iDdrhi 



d t rhi 



,y 2 dU 



(2.9) 
(2.10) 



The equations conserve the 'energy' 



Dd 
-1— mn 



U(d' 



mi H U 

3 2 1 2 



(2.11) 



The conservation law is satisfied since iCi — y 2 U/2 does not explicitly depend on t. The 'energy' 
( [2. 11| ) is equal to zero since at t — > —00 the value of rh\ should tend to zero. The property can be 
treated as the extremum condition appearing at variation of i£ r — y 2 U/2 over the initial value of <p\. 
Equating (|2~TT| ) to zero we can express m\ via 4>\. Next, since (|2.11|) is zero then the saddle-point 
value of Z(y) ( |2.8j ) can be written as i J d0i mi, where the integral over 0i goes from zero to infinity. 
Substituting into i / d0i rhi the expression of fh\ in terms of <p\ we get for the single-point statistics 



\nZ(y) 



d(d + l) 
6 



1 + 



l2y 2 P 2 



\ Dd?{d 2 -l) 



ln(LA) 



(2.12) 



Note that the expression ( |2.12|) has (as a function of y) two branch points y = ±iy sing where 

2 Dd 2 [d 2 - 1) 



/sing 



12Po 



(2.13) 



The same procedure can be done for the passive scalar difference, or, more precisely, for the object 
( p.,101 ). Taking into account the presence of the jumps ( p.7|) and ( |2.6|) in the potential U we get an 
answer slightly different from ( |2.12| ) 



\nZ(y) 



d(d + l) 
6 



2% 2 P 2 
\ x Dd\d 2 - 1) 



ln(r A) 



'Sing 



Dd 2 (d 2 - 
2APo 



(2.14) 
(2.15) 
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Note that (|2.14j) does not depend on the pumping scale L but still depend on the cutoff A. 

The characteristic value of <pi is determined by the quantity ( |2.6j ) which is much larger than unity. 
Then it follows from (|2.5| ) that exp(2pj — 2pj) < 1, i> j, (excluding for a short initial stage of the 
evolution) and we see from (|1.32j ) that fluctuations of the fields r\ are suppressed in comparison, say, 



with p a . This justifies neglecting the fields r\ and p leading to the reduced Lagrangian fl2.1|). Next, 
dynamics of <p a for a > 1 is diffusive and it follows from (|2.3|) that the characteristic value of a can 



be estimated as yDd\t\. As follows from ( |2.3| ) dt<pi ~ Dd 5 ^ 2 and we find from ( |2.6| ) the instantonic 
life time 

t lt = D^tT 2 ln(LA), (2.16) 
which determines times producing nonzero contribution to the effective action. At \t\ ~ t\ t the 



characteristic values of <p a for a > 1 are of order y\n(LA)/d and we conclude that 

4>n 1 



dJ\n(LA) 



<1, (2.17) 



at times \t\ ~ t\ t . The inequality (|2.17| ) justifies passing to the Lagrangian (|2~4T ). The same arguments 
can be applied to the generating functional for the passive scalar difference, the only modification is 
in the substitution ln(LA) — > ln(r A). 



There are also additional applicability conditions for the answers ( |2.12| , |2.14T ). To establish the 
conditions one should go beyond the main order of the saddle-point approximation. It will be more 
convenient for us to develop an alternative scheme which enables to find the conditions simpler. 
That is the subject of the next subsection. 



B. Schrodinger equation 



Here we present another way to get the answers ( [2.12| , |2T4l) . As previously we start from the path 
integral representation ( 1.33Q for the generation functional Z[y). 

Unfortunately it is impossible to get a closed equation for Z{y). To avoid the difficulty we introduce 
an auxiliary quantity 

,2 



dt' iC-^U 



(2.18) 



p(-t)=pa,ri(-t)=ri 



ty(t,y,p ,r]o) = J VpVrfDmVp exp 

It follows from the definition Q2.18] ) that 

Z{y) = lim / J] dp a d Vij V(t, y, p, 77) . (2.19) 
The equation for the function \1/ can be obtained from the expression ( |1.32|) and the definition ( p,18| ): 



d t V = 



Dd 











d 2 d d 

+2£exp(2p i -2#)— 2- +4 ^ exp(2p fc - 2 Pi ) Vkj 

oriij i<k<j oVij °Vik 



i<j 



d d 

+2 J! exp(2p fe - 2pi)- — t] km r] kn 



i<k<m,n 



drjij drji 

y 2 u T 



(2.20) 



10 



We see that the equation ( |2.20| ) for \l/ resembles the Schrodinger equation. The initial condition to 
the equation can be found directly from the definition (|2.18|) : 

*(t = 0,y,p,r ) ) = H6(p a )6(7 Hj ). (2.21) 

The value of Z in accordance with ( |2.19| ) is determined by the integral of \1/ over 7] and p. This 
integral is equal to unity at t — and then varies with increasing time t due to U ^ since only 



the term with U in ( |2.21| ) brakes conservation of the integral. Thus to find Z we should establish 
an evolution of the function \l/ from t — to large t. 

Below we concentrate on the single-point statistics. The scheme can be obviously generalized for 
the passive scalar difference. 



Let us first describe the evolution qualitatively. The initial condition (|2.21| ) shows that at t — 
the function ^ is concentrated at origin. Then it undergoes spreading in all directions except for 
pi + . . . + p d since the operator in the right-hand side of ( |2.20|) commutes with p x + . . . + p d . This is 
a consequence of the condition Det T = 1 (to be satisfied) which implies that during the evolution 
pi + . . . + pa — 0. That means that we should treat a solution of ( |2.2(J| ) which is \1/ oc 5(pi + . . . + p d ). 
The function \1/ is smeared diffusively with time and also moves as a whole in some direction, which 
is determined by the term with the first derivative in (|2.20|) . The rate of the ballistic motion is 

ln . d(d-2i + l) , 
(dtPt) = D -± >- . (2.22) 

Therefore \1/ describes a cloud, centre of which moves according to the law 

^d(d-2i + l) . 
Pi = £> — —^t. (2.23) 



Effective diffusion coefficients for rfs fall down with increasing t since in accordance with ( |2.23 ) the 



differences pk — p% [figuring in ( [2.20|) 1 are negative and grow by their absolute value. Therefore the 



diffusion over rj stops when the characteristic values of pi — pk becomes greater than unity. Note 
that the 'frozen' values of i] do not depend on y since U can be considered as uniform during the 
initial stage of the evolution. After that rj are frozen, the diffusion continues only over p's. If the 
cloud is inside the region where U ~ Pi then the evolution of the cloud is not influenced by U. After 
the period of time t\ t ( j2.16| ) the cloud reaches a barrier where the potential U falls down from Pi to 
0. The subsequent history depends on the value of y. For moderate y the cloud passes this barrier 
and continues to move with the same rate. After this the integral of \l/ will not change in time, and 
its value will determine the generating functional Z(y). Naive estimations would give the answer 
In Z(y) = —y 2 t\ t /2, which reproduces the value of the pair correlation function of 9. 

A special consideration is needed if \y\ 3> y S in g or if y is close to ±iy sing where y S i ng is introduced by 
( |2.13j ). Just the last region determined the PDF's and is consequently of a special interest. Note that 
y = ±iy sing corresponds to appearance of a bound state near the pumping boundary (where U falls 
down from P 2 to zero). If y 3> j/ S mg then the front side of the cloud reaches the jump of the potential 
much earlier than £a- The residue of the cloud (living inside the potential well) is damped due to 
the term with y and don't give a contribution to Z(y). If \y\ ^> y sing then Z(y) ^> exp(— y 2 t^/2) 1 
really the asymptotics of Z{y) is exponential in the case. If \y ± iy S i ng | y S in g then the cloud stays 
near the pumping boundary for a long time, that is the shape of \1/ inside the region U ~ P<i varies 
in time comparatively slow. Besides, a part of ^ percolates out to the region where U ~ and the 
integral of \1/ grows with increasing \t\. As y come to iy S mg closer this stage lasts longer. One can 
say that the back side of the cloud \1/ gives the right answer for Z(y). The important point is that 
if y is not very close to h/ S mg then during the time of exiting \l/ from the potential the width of \1/ 
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in terms of diffusive degrees of freedom is much less than In LA. This means that the function \]/ is 
really narrow, that justifies our consideration. 

For an quantitative analysis it is convenient to pass to the variables fa ( |2.2| ). Since an //-dependence 
of \I/ is frozen after an initial part of evolution, then it is possible to obtain an equation for the integral 
of \& over rj: 



where we included also an integration over 
for the function (|2.24j ) is 



d0 rf n dr fe^ ( 2 - 24 ) 
b d to remove the factor 8(p\ + . . . + pa). The equation 



Dd 
IT 



d-l Q2 



ld(d 2 -l) d 



r y 2 u~ T 



(2.25) 



where U is function of <p a only which can be found by substituting into U the 'frozen' values of 77's. 
Qualitatively U has the same structure as U itself. One can conclude from ( |2.25f) that the cloud 
described by \& moves ballistically into the fa direction and spreads along other directions. We are 
going to treat the situation when the cloud remains narrow during the relevant part of the evolution. 
Then one can integrate \I/ over all fa, i > 1 in a similar way as in the case with 77's and get an Id 
equation for 



i-l 



Tldfa*. 



The function \1/ satisfies the equation 



a t Dd 

dt^ = — 



d 



ld{d 2 - 1) 



d - y 2 U - 
— — - — \1> . 



(2.26) 



The initial condition to Eq. ( |2.26| ) is ^f(t — 0) = 5(<pi). The potential U is obtained from U by 
substitution <p a — > for a > 0. In fact, on the direction ( |2.23| ) we have that is the potential U depends 
only on p\. The barrier is achieved when p% ~ In LA. Passing to the variables fa we conclude that 
the potential U diminishes from P 2 at fa < 0a to zero at fa > 4>a where 0a is defined by Q2.6Q . 

The character of a solution of the equation ( |2.26| ) can be analyzed semiqualitatively in terms of 
the width I of \I/ over fa and its amplitude h. When ^ reaches the pumping boundary then it stops 
there for a period of time. Then the width I and the amplitude h governed by the equations 



d/ r,^ Dd 

dt = - DdX+ - 



dh 
dt 



Ddh y 2 P 2 h 



I 2 



(2.27) 



where A = yd(d 2 — 1)/12, DdX is the rate of the cloud motion along fa direction (when U = const), 
and Dd is the diffusion coefficient for fa direction. One can estimate from the 1st equation the width 
I as I ~ 1/A. Then from the 2nd equation the height h falls down or grows in time depending on y. 
The characteristic y where the regime changes is of the order 1 2/ S ing | 2 ~ DdX 2 / P2. We'll show this by 
consistent calculations. 

The equation ( 2.26|) can be solved analytically, e.g., by the Laplace transform over time t. Making 
the Laplace transform one gets 



ptf (p) - 5(fa) = ^ 



_d_ 



'd{d 2 - 1) 



t (p) - y-u(fa)^(p) 



(2.28) 
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We are interesting in the bound state describing by the equation. Solutions for *f?(p) in the intervals 
(—00, 0), (0, 0a), (0 a, 00) are exponential, one should match them. The function ^f(p) as a function 
of p has two branch points at 

Dd 2 (d 2 - 1) - . Dd 2 (d 2 - 1) 
Pi = YA ~~ y Ps/2 ' P2 = 24 ' (2 ' 29) 

coming from the regions 0i < a and 0i > 0a, respectively. When one of this branch points passes 
p = then ^ starts to grow exponentially in time. This happens when y passes ±iy sing from fl2.13|) 
moving along the imaginary axis. 

The value of the generating functional is determined in accordance with fl2.19|) by large time 
behavior of ^(t). This means that we should be interested in the behavior of ^(p) at small p. Really 
the function / d0i ^(p) entering ( |2.19| ) has a pole at p = related to the asymptotic behavior 



( 2p 
^(p) oc exp — 



Dd\] d(d 2 - iy L j ' 

at 0i > 0a and small p, the behavior can be found from ( [2.28|) . Just the residue of / d0i ^f(p) at 
the pole determines Z(y). To find the value of the residue we should analyze the behavior of ^>(p) 
at < 0i < 0a- At small p there are two contributions into ^ which behave as 



-exp .p^ ± ,^_±Z + iL^U , (2.30) 




as follows from (|2.28| ) at p — 0. Therefore the residue which is determined by the integral / d0i ^(p) 



over the region 0i > 0a is proportional to 




Substituting here (|2.6|) we reproduce (|2.12j ). 



Let us now establish the applicability condition of the above procedure. The expression ( [2.31|) 
implies that the exponent with the sign minus in ( |2.30| ) gives a negligible contribution to at 
0i = 0a- The condition is satisfied if 



\y 2 + yL g \<Pl^ir 
±2 



Substituting here (|2.6| , |2.13|) we obtain 



smg 



/sing 



> (d 4 ln 2 LA) 1 . (2.32) 



For y closer to ±iy S i ng one should be careful since then a fine analytical structure of Z(y) will be 
relevant. As an analysis for d = 2 shows |25j Z(y) has a system of poles along the imaginary semiaxis 

starting from ±k/ S i ng and the parameter (d 4 \n 2 LAj just determines the separation between the 
poles. The assertion about the cut made in the previous subsection is related to the restrictions of 
the saddle-point approximation which cannot feel this fine pole structure, it gives a picture averaged 
over the interpole distances, it is precisely the cut. This approach is correct just at the condition 
(PI). 
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Note that the same criterion ( |2.32| ) justifies our assumption that the cloud described by $ is narrow 
during the relevant part of the evolution. Namely, the duration of the part is determined by the 
time tgxit — Pi 1 (see ( 2.29| )). This is the time which the cloud stays near the barrier. For y close to 
±^/sing the time can be estimated as t~^ t ~ -P2 1 2/sin g | \y T^singl- Then the diffusive width \/Ddt exit 
of \I/ in the directions <j) a for a > 1 is much less than <f)\ just if ( |2.32| ) is satisfied. Principally the 
diffusive dynamics at d > 2 could modify the noted fine pole structure of Z, the problem needs a 
separate investigation. 

The same procedure can be done for the passive scalar differences. The cloud \1/ should pass 
the region p\ < ln(L/ro) before it reaches the potential. Then it enters the region U = 2P2 with 
some finite diffusive width. One can note, however, that it is irrelevant. The only characteristics of 
potential which are needed are its value (here 2P2 instead of P2) and the length of the path inside it 
(which is Api = ln(roA) instead of ln(LA)).The evolution of \I/ goes in the same way as in the case 
of single-point statistics. Again, we get the answer 02.14j ) and the criterion analogous to (|2.32j ). 

In the subsection we presented the analysis based on the dynamical equation (|2.20| ) for the auxiliary 
object The results obtained can be reproduced also on an alternative language: For this we should 
introduce another auxiliary object the equation for which is stationary. The corresponding scheme 
which could be interesting from the methodical point of view is sketched in Appendix. 



III. CALCULATION OF PDF 

In this section we calculate the PDF's V for the objects (|1.8|) and (|1.10|) . The most convenient 
way to do it is in using the relation 

I p- eM-iy#)z(y) , (3.1) 

J 27T 

where $ is 

<&= J dr/?(r)0(O,r). (3.2) 
Let us remind that knowing V($) one can restore also moments of 

l )= [ dtf |tf| n 7>(tf). (3-3) 



The generating functional in (|3.1| ) is determined by ( [2.12j ) or ( 2.14J) . Being interested in the 



main exponential dependence of the PDF's for the objects (|1.8|) and (|1.10| ) we can forget about 
preexponents. Then 

dy 



P(0) = I ^exp {-jy# + q [l- y /l + y^/y^ ) , (3.4) 



where 



Dd 2 (d 2 -l) /r . 2 Dd 2 (d 2 -1) 
V S m R = j^p ( for scalar), y sing = — (for difference), (3.5) 



q = ^ ^~ — - ln(LA) (for scalar), q = — - ln(roA) (for difference). (3.6) 



Since both q defined by ( |3.6| ) are regarded to be much larger than unity the integral ( |3.4| ) can be 
calculated in the saddle-point approximation. The saddle-point value is 
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Then 

/ 

lnP(i?)~g 1- 
V 



\ 



1 + I • (3-8) 



The expression leads to the exponential tail 

lnP(tf)~-y sing |tf|, (3.9) 

realized at \$\ ^> q/y S ing- The coefficient y S i ng in ( |3.9| ) determined by ( |2.13| ) is in accordance with 
the result obtained in [|i~2]j . 
The expression ( |3.8|) enables one to find the following averages in accordance with ( |3.3j ) 

2 P AP 

^ = w^ HLK) - {m ' )2) = w^m Hr '' Kh (310) 



The expressions ( j3.10|) can be obtained also by direct expansion of Z(y) from (|2.12|) or ( 2.14 ). The 
universal tail (13.91) is realized if 



9 A » dln(LA) , A# A » V((A0 A ) 2 ) dln(r A) . (3.11) 

Since both logarithms are supposed to be large we conclude that there exists a relatively wide region 
where the statistics of d is approximately Gaussian, the region is determined by the inequalities 
inverse to (ft.llj) . 

Let us discuss the applicability conditions of the expression ( |3.8| ). First, if one calculates the 
passive scalar PDF by the saddle point method then the position of the saddle point is determined 
by ( ^32|) if 



l 9<d 2 y^ln 2 (LA). (3.12) 

The applicability domain of the saddle-point method overlaps the region of validity of expression 
fl2.12j ) for the generation function Z(y). The above inequalities are correct for 8\, for A8\ one should 



substitute ln(LA) by ln(r A). Second, fluctuations of y have to be small compared to the distance 
between y sp and y S mg- This gives the same criterion ( |3.12|) 



Let us stress that though formally our procedure is incorrect at $ ^ d 2 y P 2 /D\n 2 (LA) the answer 
will be the same: the PDF will be determined by the exponential tail ( |3.9| ). The point is that the 
character of the integral ( |3.1| ) at such extremely large $ will be determined by the position of the 
singular point of Z(y) nearest to the real axis. This is just h/ S i ng leading to (|3.9p. To conclude, only 



the character of the preexponent in V{$) is changed at $ ~ d 2 \J P 2 j D\n 2 '(LA) whereas the principal 
exponential behavior of V($) remains unchanged there. 



CONCLUSION 



The single-point statistics of the passive scalar 9 and the statistics of its difference A# are tradi- 
tional objects which carry an essential information about correlation functions of the passive scalar 
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in the convective interval. We examined the passive scalar in the large-scale turbulent flow where the 
correlation functions logarithmically depend on scale. Since really the logarithms are not very large 
it is useful to have the whole PDF's of 9 and AO. That was the main purpose of our investigation 
which was performed in the frame of the Kraichnan model. The single-point PDF for the passive 
scalar and the PDF for the passive scalar differences can be obtained from ( |3.8| ) if to substitute 
A — > r^jf where r<iif is the diffusive length. Though both the advecting velocity and the pumping 
force in the Kraichnan model are considered as 5-correlated in time we hope that our answer are 
universal that is are true in the limit when the size of the convective interval tends to infinity for 
arbitrary temporal behavior of the velocity and of the pumping. The reason is that the spectral 
transfer time grows with increasing the convective interval and in the limit is much larger than the 
correlation times of the velocity and of the pumping. 

We believe also that the analytical scheme proposed in our work could be extended for other prob- 
lems related to the passive scalar statistics. Note as an example the work [T3| where a modification 
of the scheme enabled to find the statistics of the passive scalar dissipation. It is also useful for in- 
vestigating the large-scale statistics (on scales larger that the pumping length) of the passive scalar 
p6fl . We hope also that it is possible to go beyond the case of the large-scale velocity field using the 
perturbation technique of the type proposed in |f27j-|29ll. 
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APPENDIX A: 



Here we present an alternative way to obtain the answers ( |2.12| , |2TT41 ). We will use an auxiliary 
quantity 



Z(y,po,Vo) = J VpVrfDmVp exp J dt (i£ - yf7 



p(o)=po,v(o)=vo 



then 



Z(y)=E(y,0,0). 

Function H can be also defined as 

Po, Vq) = lim / n dp a dr^ tf(t, y, p, rj) 



(Al) 



(A2) 



(A3) 



where \I> is governed by equation ( |2.20|) with initial condition ty(t = 0, y, p, rj) = S(p — Po)5{rj — rj ). 
The equation for the function H can be found from ( |1 . 3 2| , [AT]) : 



d Q2 



d 



d 



t d Pi d \7^i d P 



i=l 
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d 2 3d 
+2£exp(2p j -2 Pi )— ¥ + 4 £ exp(2p fc - 2 Pi )r) kj — — 
i<j dr lij i<k<j oriijdriik 



d d 

+2 2^ exp(2p fc - 2pi)r\ km r\ kn - — 



The boundary condition for the equation (|A4j) follows from the definition (|AT|) : For large enough 
Pi, rji the potential U = for i = and remains zero also at finite times t. Therefore the integral 
( |AT| ) should be equal to unity in the case. Thus E(y,p,rj) should tend to unity where p, rj — > oo. 
Let us rewrite the equation (A"4) in terms of the variables (|2.2|) : 



(fx + 7) (Si + Q = 0, E = Si + £ , (A5) 
~ <9 2 /d(d 2 - 1) d y 2 U 

r i = 7^2+1 — q — ^T-^T' ( A6 ) 



d-l ^2 ^2 9 9 

7 = H + 2^ ex P( 2 Pfc - 2 Pi)^~T~ + 4 exp(2p fc - 2pi)rj kn 



,2'"^ , .^ <n <),,,„ <),,,, 

_d d_ 

, i, ,„_,, dr] im drji, 



d d 

+2 ^ exp(2p fc - 2pi)rj km rj kn - — . (A7) 



Here U as a function of <fii is equal to P2 inside a region restricted by </>X and <p\ (where 0^ are 
functions of variables <p 2 , ■ ■ ■ , <frd, if) and tends to zero outside the region. We will solve the equation 
using the perturbation theory over 7, £. Then the zero order equation is 



TxSx = . (A8) 
The equation ( |A8|) can be easily solved at 0^ < 0x < 0^, the answer is 

2A 



A 2 + ^ + A 




exp - \ M 2 + ^-A (0j-0i) , (A9) 



where A = Jd(d 2 — 1)/12, DdX is the rate of the cloud motion along 0i direction. The result flAU 



can be obtained using the inequality J X 2 + y 2 P2/ Ddln LA 1. The derivative dEi/d^i = at 
0i < 4>\- However, Si ^ 1 in this region. This is due to the following fact: this region corresponds 



to the evolution of \I> when its initial position is to the left of potential U (see ( |A3f) ). During the 
evolution cloud ^ passes the region of U and its integral over p, 77 changes. Then S is not equal to 
1. Only when the distance between initial position and potential is of order ln 2 LA the diffusion of 
the cloud leads to smallness of the part of \I> which passes the potential U, and H becomes closer 
to unity. Thus, function H has long tail from potential with direction to negative 0i where it is 
not equal to 1. The procedure of finding 5 from the equation (|A8|) corresponds to the geometrical 
optics approximation (taking into account only derivatives in propagating direction; this allows one 
to get the fact of propagation). This tail of S in this approximation is nothing else than the shadow 
of potential U. Higher orders of perturbation theory over the transverse derivatives correspond to 
diffraction corrections. 

Now let us consider the correction £. The equation for it looks like (r\ + 7)^ = — 7S1. Again let 
us neglect 7 in l.h.s. and solve the equation. Si is some exponential function with scale of the order 
1. Then 7S1 ~ Si. Not that 7E1 is almost equal to zero at 0i > To estimate £ one should 
construct Green function (L7(0i|0 o ) for operator T\\ 
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- I - I \I X2 + ^TTT - A Uo I I 1 - Oxp ! -2WA 2 + ^§(0+ - 0o) I I , (A10) 



2A \ IV Dd I I \ IV Dd 

where 




C = (y^ 2 + y 2 P 2 /Dd - a) / (^A 2 + y 2 P 2 /Dd + A 



The unity in brackets in ( |A10| ) gives the correction for S which has the same exponential factor 
as Si. Thus £ does not change the answer with the logarithmic accuracy. The second term in the 
brackets gains while 0o is close to cf)\. It is due to nonzero width of the cloud 5 and to dependence 
of tit on other variables. Again, it does not change the exponent. 
The case of the passive scalar differences can be considered in a similar way. 
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